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A class of braneworld models can encounter a 'quiescent' future singularity 



Braneworld models of dark energy are examined in the light of observa- 
tions of high redshift type la supernovae. Braneworld models admit several 
novel and even exotic possibilities which include: (i) The effective equation 
CO \ of state of dark energy can be more negative than in LCDM: w < —1] (ii) 

o , 

psj , at which the energy density and the Hubble parameter remain well behaved, 

o ^ 

fH \ but higher derivatives of the expansion factor (a, a etc.) diverge when the 

^: 

Q , future singularity is reached; (iii) The current acceleration of the universe is a 

H , 

' transient feature in a class of models in which dark energy 'disappears' after a 

certain time, giving rise to a matter dominated universe in the future. Since 
^ ' horizons are absent in such a space-time, a braneworld model with transient 

acceleration might help reconcile current supernova-based observations of an 
accelerating universe with the demands of string/M-theory. A maximum like- 
lihood analysis reveals that braneworld models satisfy the stringent demands 
imposed by high redshift supernovae and a large region in parameter space 
agrees marginally better with current observations than LCDM. For instance, 
models with w < —\ {> —1) provide better agreement with data than LCDM 
for ^Irr, > 0.3 ( < 0.25). 
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I. INTRODUCTION 



The idea that the universe may have more than three spatial dimensions first appeared in 
the works of Kaluza and Klein over 70 years ago. The notion of extra dimensions gradually 
became popular because of the hope that one might succeed in relating the gauge symmetries 
of particle physics to the isometries of a compact higher dimensional manifold [1] . All such 
models, however, assumed that the compactification scale is small (/ ~ Ipi ~ 10~^^ cm) 
hence unobservable. More recently, the many variants of superstring theory allow for the 
possibility that at least some of the extra dimensions of nature may be macroscopic. For 
example, the eleven -dimensional supergravity model of Horava and Witten [2] assumes that 
ordinary matter fields are confined to a submanifold (brane) which is embedded in a higher 
dimensional space (bulk). In an important recent development, Randall and Sundrum (RS) 
examined a simplified variant of this model consisting of a three dimensional brane embedded 
in a four dimensional anti-de Sitter (AdS) bulk [3] . Their results showed that gravitational 
excitations are confined close to the brane giving rise to the familiar law of gravity, and 
suggesting that one could identify the brane as our observable universe. 

Subsequently the RS ansatz was generalised to incorporate both expanding FRW-type 
models [4] and anisotropic space-times [5]. Of particular importance to the present study 
is the observation that simple extensions of the RS scenario can give rise to a universe 
which is accelerating, in agreement with studies of high redshift supernovae [6,7]. It is now 
well established that high redshift type la supernovae appear fainter than expected in a 
spatially fiat matter dominated (Einstein-de Sitter) universe [8,9]. One way of explaining 
this discrepancy is to postulate that the universe is filled with a smooth component carrying 
large negative pressure (dark energy). Although several possible candidates for dark energy 
have been suggested (the cosmo logical constant, quintessence etc.) none is entirely problem 
free (see [10,11] for recent reviews). 

In this paper we shall focus on a new form of dark energy based on the braneworld model 
examined in [7,12] (see also [14,15]). Braneworld models of dark energy have interesting new 
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properties including the fact that, depending upon the form of bulk- brane embedding, the 
effective equation of state of dark energy can be i/; > — 1 or w < — 1. In addition, for an ap- 
propriate parameter choice, the acceleration of the universe can be a transient phenomenon, 
thus helping reconcile high-z supernova observations of an accelerating universe with the 
requirements of string/M-theory. 



II. DARK ENERGY FROM BRANEWORLD MODELS 



The equations of motion governing the braneworld can be derived from the action [14,15] 



[ (i?5 - 2Ab) - 2 / k\+ [ (m^Ri - 2^) + / L (V, 0) ■ (1) 

Vbulk ./brane J J brane ^ brane 



Here, is the scalar curvature of the metric Qah in the five-dimensional bulk, and R^ is the 
scalar curvature of the induced metric h^a on the brane. The quantity K = Kaph""^ is the 
trace of the extrinsic curvature on the brane defined with respect to its inner normal. 
L{hai3, (t>) is the four-dimensional matter field Lagrangian, M and m denote, respectively, the 
five-dimensional and four-dimcnsional Planck masses, Ab is the bulk cosmological constant, 
and (7 is the brane tension. Integrations in (1) are performed with respect to the natural 
volume elements on the bulk and brane. 

The presence of the brane curvature term /brane -^4 (1) introduces an important 
length scale into the problem I — 2m? /M^. On short length scales r <^ I (early times) one 
recovers general relativity, whereas on large length scales r ^ / (late times) brane-specific 
effects begin to play an important role, leading to the acceleration of the universe at late 
times [6,16,7]). 

The cosmological evolution of the braneworld is described by the Hubble parameter 



^2 ^ _^ ^ p + a ^ 2_ 
a? 2>'m? P 



1±\ l+l 



p + a _ ^ _ £ 
3m? 6 



(2) 



The two signs in (2) correspond to the two separate ways in which the brane can be embed- 
ded in the higher dimensional bulk. The underlined term in (2) makes braneworld models 



different from standard FRW cosmology. The limiting case of our model m — corresponds 
to the well known FRW generalisation of the RS scenario [4] 



In this case, braneworld evolution departs from the standard FRW law at early times when 
p/a 3> 1. However, as remarked earlier, braneworld models described by (2) depart from 
FRW behaviour at late times, a property that opens radically new possibilities for the present 
and future state of our universe. 

Braneworld models fall into three main categories [7]: 

• 1. BRANEl (Bl): The Hubble parameter is given by 

-"o 

(4) 



where z — ao/a{t) — 1 is the cosmological redshift, while 
o Po o - — o ^ o - ^— o 

are dimensionless parameters whose values must be determined from observations (the 
subscript '0' refers to their current value), fl^^ is determined by the constraint relation 



n^ + n^ + n^- 2^ni + = i. (6) 

2. BRANE2 (B2): The Hubble parameter is given by 

^Qfi = n^{i+zf + n,{i+zf + + 2^1 + 2^iJn^{i+z)^ + + + Qa, , 

(7) 

where fli < 1 + fl\^ and fl^^ is determined from 



n^ + n^ + n^ + 2JniJi-n^ + nA^ = i. (8) 



The two models BRANEl and BRANE2 are complementary and reflect the two dis- 
tinct ways in which the brane can be embedded in the bulk [7]. Clearly, by setting 
fli = the underlined terms in (4), (6), (7), (8) vanish and we recover the LCDM 
model in both cases. The underhned terms in (4) & (7) are caused by braneworld 
effects and give rise to an important general result, namely 



H{z) 



Bl 



<H{z) 



, H{z) 

LCDM 



<H{z) 

LCDM 



B2 



<H{z) 



(9) 

SCDM 



ie the universe expands at a faster (slower) rate than LCDM in BRANE2 (BRANEl). 
Since most cosmological observables involve H{z) either directly or indirectly, 
braneworld models can exhibit properties which can be quite distinct from those of 
either LCDM or SCDM. 

• 3. Disappearing dark energy (DDE): The Hubble parameter is given by (7), and 
the cosmological parameters flm, ^a,^i, are constrained to satisfy the following 
relations [7]: 

< Qa, ^ Qa, > , (10) 



0^ + 2^(^1 + 0a, -y^) = 1, (11) 



Q,^-2^ninA^. (12) 

(We make the assumption that the universe is spatially flat, so that fl^ = 0.) 

From (4) & (7) it is easy to see that all braneworld models approach the standard matter 
dominated universe at early times [with a small correction term ~ (1 + z)^/^]. At late times 
the behaviour of the braneworld can differ from both LCDM and SCDM. This feature makes 
braneworld models testable and allows the braneworld scenario to provide a new explanation 
for the observational discovery of dark energy. 
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The braneworld models described above provide a common platform for understanding 
the properties of dark energy. For instance the expression for the Hubble parameter in these 
models (4) & (7) allows us to explicitly determine key cosmological quantities including: 



• The luminosity distance dL{z): 



dL{z) 
1 + z 



Jo 



dz' 



(13) 



the angular-size distance 



dA{z) 



dz' 



l + zJo H{z') ' 



(14) 



• the deceleration parameter: 



q{z) = 



H'jz) 
H{z) 



(15) 



the effective equation of state of dark energy: 

2q{z) - 1 



w{z) 



3[i-n^{z)] ' 



(16) 



• the age of the universe: 



"^)=/ (TT7 



dz' 



(1 + z')H{z') ' 



(17) 



the 'statefinder pair' [17]: 



= 1 



am 
r — 1 

3(g-l/2)' 



H" fH'' 



[l + zf-2^{l + z), 



(18) 



the product dA{z)H{z), which is used in the Alcock-Paczynski anisotropy test [18], 



the product d\{z)H ^{z), which plays a key role in the volume- redshift test [19]. 



While providing a common basis for the existence of dark energy, braneworld models 
Bl, B2, DDE, have important properties and attributes which serve to distinguish these 
models both from each other, and from alternate models of dark energy such as LCDM 
& quintessence. For instance, the luminosity distance in BRANEl can be larger than the 
luminosity distance in LCDM: d^lBi > di\LCDM- This follows from (9) and leads to an 
important result. Namely, using (16) we find the current value of the effective equation of 
state of dark energy: 



from where we see that wq < —1. BRANEl therefore has the remarkable property that the 
effective equation of state of dark energy can be more negative than that associated with a 
cosmological constant {w — —1). This feature distinguishes Bl models from LCDM as well 
as from scalar field based quintessence models. 

BRANE2 has the opposite property, namely dilscDM < c^l|s2 < dLlicDM, which trans- 
lates into 



We therefore see that Bl & B2 are complementary models and that the effective equation 
of state in Bl (B2) is softer (stiffer) than the w — —1 typical of LCDM. Furthermore in 
both B2 & DDE, the current acceleration of the universe - like earlier matter and radiation 
dominated epochs - can be a transient feature. In the case of DDE, the current accelerating 
phase will be replaced by a matter dominated epoch during which ty ~ [7] . In Fig. 1 we 
show the behaviour of the deceleration parameter for this class of models. It is seen that 
braneworld dark energy disappears in the future and the universe becomes matter dominated 
as z — > — 1. It is easy to show that the present value of the effective equation of state of 
dark energy in DDE is given by 




(19) 



-1 < Wo < 0, 



where 




(20) 




(21) 
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FIG. 1. The behaviour of the deceleration parameter with redshift for DDE. Sohd curves from 
top to bottom are DDE models with = 0.2, and ^Ia^ = 1.0, 1.5,2.0 respectively. The dashed 
line is SCDM with $7^ = 1, the dot-dashed curve is LCDM with = 0.2. The vertical dotted line 
marks the present epoch at z = 0, and the horizontal dotted line marks a q = Milne universe. In 
the DDE models, the universe ceases to accelerate and becomes matter-dominated in the future, 
unlike the LCDM model which remains dark energy dominated at all future times. 

Since Qi < Qa^ we find that wq > —0.5. Model B2 permits more exotic possibifities. In a 
subclass of models acceleration gives way to an epoch during which the universe decelerates 
at an increasingly rapid rate [12]. The expansion of the universe culminates in a 'quiescent 
singularity' which is distinguished from conventional general relativistic singularities by the 
fact that the energy density and Hubble parameter remain finite, while higher derivatives of 
the scale factor (a, a etc.) diverge, when the 'future singularity' is reached (Fig. 2). (The 
future singularity, measured from the present epoch, is reached in a few billion years for 
most B2 type braneworld models [12].) 

Since neither B2 nor DDE possess an event horizon both can successfully reconcile the 
demands of string/M-theory with a universe which is currently accelerating [13]. 
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FIG. 2. The BRANE2 universe can encounter a singularity lying in the future as demonstrated 
in this figure. The deceleration parameter becomes infinite as the singularity is approached while 
the Hubble parameter (and the density, pressure) remain finite. The vertical dot-dashed line 
corresponds to the present epoch 2 = while the dashed line represents the dimensionless Hubble 
parameter. The solid line shows the deceleration parameter q{z). The model parameters are 
^rn = 0.15, J7; = 0.4. Although permitted by supernovae observations this particular model 
appears to be disfavored by clustering bounds on Qrn (see Fig. 6). 
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III. METHODOLOGY AND RESULTS 



We constrain the parameter space of braneworld cosmology by ensuring that our cosmo- 
logical models provide a good fit to Type la supernova data. For this purpose we use the 
54 SNe la from the primary 'fit C of the Supernova Cosmology Project, which includes 16 
low redshift Calan-Tololo SNe [9] . Fit C is a subsample of a total of 60 SNe of which six 
are excluded as outhers: two low redshift SNe due to suspected reddening and four high 
redshift SNe of which two are excluded due to atypical light curves and two because of 
reddening. The measured quantity in this data, the bolometric magnitude m^, is related to 
the luminosity distance and therefore the cosmological parameters by the following equation 



where Dl = HqcIl is the Hubble-parameter-free luminosity distance and M. = Mb + 25 — 
5 logio-f^o is the Hubble-parameter-free absolute magnitude. We shall assume that the SNe 
measurements come with uncorrelated Gaussian errors in which case the likelihood function 
is given by the chi-squared distribution with N — n degrees of freedom C oc exp (— 
( In our case N — 54, and n — A: Bl and B2 estimate four parameters each, and DDE 
estimates three parameters with one constraint, Eq (10) lowering the degrees of freedom) 
The x^-statistic is defined as 



where mf is the effective B-band magnitude of the i-th supernova obtained after correcting 
the observed magnitude at redshift z for the supernova width-luminosity relation, (7^. is the 
error in magnitude at redshift z, and m{zi) is the apparent magnitude of the i-th supernova 
in the braneworld model. 

For BRANEl and BRANE2, the parameters to be estimated are M., 0,^, ^i, {^a 
is calculated from Eqs (6) and (8) respectively for Bl & B2). For our purposes the quan- 
tity is a statistical nuisance parameter, and we marginalise over it assuming no prior 



niB = M + 5 logioDiiz; Q^n, ^i, ^^aJ , 



(22) 




(23) 
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m 

FIG. 3. The likelihood function is shown as a function of each of the parameters ^1^, ^i, i^Aj, 
and M after the remaining parameters have been marginalised. The top, middle and bottom panels 
show the likelihood curves for Bl, B2, and DDE. The solid lines correspond to the value of the 
likelihood function for a given parameter after it has been marginalised over all other parameters, 
while the dashed lines show the likelihood function evaluated by fixing the other parameters at 
their maximum likelihood value. In all three cases the likelihood function is normalised with the 
maximum likelihood value set to unity. 
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knowledge to get the three-dimensional probability distribution in the (Q^,Q;,QaJ space: 
P{^m, ^aJ = / -P(^m, ^Afc) ■M)dM.. We perform maximum hkelihood analysis on the 
system with the priors < flm < 1, >0, fl\^ > 0. For B2, we use the added constraint 
Qi < 1 + Qa^. For DDE, the parameters to be estimated are M., ilm, ^a^- (The remaining 
two parameters and fla- are related to flrn, through Eqs (11) and (12) respectively). 
For this model we use the prior < flrn < 1 and the constraint Eq (10) (see [7] Appendix). 
For Bl & B2 the constraint relations (6), (8) combined with Q„ = 0, set the lower bound 
> — 1- However since fl\^ > is a more physically appealing model (it includes anti-de 
Sitter space (AdS) bulk geometry), we choose this as a prior for further analysis. Results 
for — 1 < Qaj < will be presented elsewhere. 

In Fig. 3 we show the hkelihood curve as a function of Qm, Qi, f^A^ and M. for the 
braneworld models Bl, B2, and DDE. We see that for all three cases the likelihood is a 
sharply peaked function of A4, which is non-zero over a very limited range. Therefore it 
seems reasonable to marginalise over A4 in this range. Surprisingly, the likelihood function 
for BRANEl peaks at f^Aj = 0, which is its best-fit value. Thus the BRANEl universe 
appears to prefer a vanishing cosmological constant in the bulk. The BRANE2 universe 
shows very different behaviour. In this case the hkelihood function is fiat and therefore 
insensitive to the value of VL^^ over the wide range < VL\^ < 5. Beyond JIa;, ~ 5 the 
likelihood function drops steeply. In Figs 4 and 6 we have shown maximum likelihood 
contours for Bl and B2 obtained by (a) marginalising over f^Ai,, and (b) setting QAt = 0. The 
close similarity between these two figures suggests that the likelihood contours marginalised 
over ilAfe do not change the results appreciably, either for Bl or for B2. 

In Fig. 4 we show confidence levels in the Q„i — Vti plane for BRANEl. We find that 
the BRANEl model can be definitely excluded only if the matter density is small Vt^ <^ 0.2. 
For Q,Tn ^ 0.3 the BRANEl model agrees well with supernovae data, with the agreement 
extending to larger values of Qi as Qm increases. Clearly in order to be able to restrict 
the braneworld parameters further one needs additional information about the dark matter 
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FIG. 4. Confidence levels at 68.3% (light grey inner contour) 95.4% (medium grey contour) and 
99.73% (dark grey outer contour) are shown in the J7;-^m plane for BRANEl. Panel (a) represents 
confidence levels in the O^-Q^ plane when marginalised over Qa;, and Ai, while panel (b) shows 
the confidence levels marginalised over A^, with = 0, the best-fit value. We see that taking the 
best-fit value of = instead of marginalising over it does not change the results appreciably. 
The dotted region represents the intersection of the 3a confidence level with the observational 
constraint 0.2 < < 0.5. The thick solid line in (b) represents the most likely value of if the 
value of Cljn is known exactly. We see that the BRANEl model is in good agreement with SNe 
observations if the value of Q^n is moderately high: ^ 0.3. (It should be noted that the value 
of the five dimensional Planck mass corresponding to $1/ 1 is M 100 MeV.) 
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(a) (b) 

FIG. 5. The shaded region represents the intersection of the 3a confidence region in the i}rn~^i 
plane for Bl with the observational constraint 0.2 < 0,^ < 0.5 (J7a(, = is assumed). In panel (a) 
we show the values of different cosmological quantities determined at the present epoch. The dashed 
lines correspond to the current effective equation of state of braneworld dark energy for BRANEl 
models: wq = —1.70, —1.50, —1.30, —1.15, and — 1.05 (top to bottom). The dotted lines corre- 
spond to the current age of the BRANEl universe: iJo*o = 0.85, 0.90, 0.95, 1.02, and 1.10 (top to 
bottom). This corresponds to io(Gyrs) = 11.8, 12.6, 13.3, 14.3 and 15.4 if Hq = 70 km/sec/Mpc. 
In panel (b) the dashed lines correspond to the epoch, za, at which the braneworld universe (Bl) 
begins to accelerate: za = 0.45, 0.60, 0.75, 0.90, and 1.05 (top to bottom). 
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density not included in the supernova observations. Currently there is no firm consensus on 
the value of Q^- While recent studies of galaxy clustering indicate Qrn ~ 0.3 [20], larger 
values of flm may be favoured by observations of clusters of galaxies [21,22]. In this paper we 
shall assume the weak clustering bound 0.2 <, fi^ <, 0.5 and study the braneworld models in 
greater detail in this region. In figure 4 the intersection of the 'Scr' SNe bound and the bound 
on flrn is shown as a dotted region. Interestingly a large region in BRANEl parameter space 
is seen to satisfy both the supernova constraints as well as the clustering bounds on Q^. 

Fig. 5 highlights the ilm — region permitted both by the SNe bound on BRANEl and 
by the clustering bound 0.2 <, flrn ^ 0.5 for the best-fit value — 0. The lines running 
through the region depict, respectively, (i) different values of the current effective equation 
of state. Wo, of dark energy, (ii) the age of the universe, to, (hi) the acceleration epoch, za, 
defined as the redshift at which the universe began to accelerate. We find that the age of the 
universe is constrained to lie in the range 0.85 <, toHo < 1.10, while the effective equation 
of state is -1.70 <, wo <, - 1 (One should note that wq < -1 in the BRANEl model) . 
From the figure we find that the acceleration of the BRANEl universe is a fairly recent 
phenomenon, which commenced at 0.45 za < 1.05. 

Fig. 6 shows the confidence levels in the Qm — ^/ plane for the BRANE2 model. Again 
the intersection of the '3cr' SNe bound with the bound 0.2 < fl^n < 0.5 is shown as the 
dotted region. In contrast to BRANEl lower values of Jl^ agree better with observations 
for BRANE2. The region to the right of the dashed fine in the figure corresponds to the 
class of braneworld models having future singularities. It is interesting to note that such 
models are in good agreement with SNe data for low values of Q^, but are ruled out (at 3a) 
if ilm > 0.25. 

Fig. 7 shows lines of constant wo, to, and za within the allowed region for BRANE2. 
Our results show that, the age of the universe, its effective equation of state, and the 
commencement of acceleration epoch, are constrained to lie within the following intervals: 
0.88 <, Hoto < 1.07, -1 <, Wo <, - 0.6, and 0.4 < za < I- (One should note that Wq < -1 
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FIG. 6. Confidence levels at 68.3% (light grey inner contour) 95.4% (medium grey contour) 
and 99.73% (dark grey outer contour) are shown in the Jlj-^^m plane for the BRANE2 model. The 
figure (a) represents confidence levels in the 0^-0^ plane when marginalised over and Ai, 
while the figure (b) shows the confidence levels marginalised over M, with CIai = 0. We see that 
fixing the value of f^Aj to zero instead of marginalising over it has negligible effect on the results. 
The dotted region represents the intersection of the 3a confidence level with the observational 
constraint 0.2 < < 0.5. The thick solid line in (b) represents the most likely value of if 
the value of is known exactly. The region to the right of the dotted line in (b) represents B2 
universes which will encounter future 'quiescent' singularities. We see that the BRANE2 model is 
in good agreement with SNe observations for lower values of < 0.30. 
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(a) (b) 

FIG. 7. The shaded region represents the intersection of the 3a confidence level of BRANE2 
for = with the observational constraint 0.2 < < 0.5. In (a) we show the values 
of different cosmological quantities determined at the present epoch. The dashed lines corre- 
spond to the current effective equation of state of braneworld dark energy for BRANE2 models: 
Wo = -0.62, -0.75, -0.85, -0.92 and - 0.98 (top to bottom). The dotted lines correspond to 
the current age of the BRANE2 universe: Hoto = 0.88, 0.92, 0.97, 1.02, and 1.07 from (top to 
bottom). This corresponds to to (Gyrs) = 12.3, 12.9, 13.5, 14.3, and 14.9 if ifo = 70 km/sec/Mpc. 
In (b) the dashed lines correspond to za = 0.40, 0.55, 0.70, 0.85, and 1.00 (top to bottom), where 
ZA is the epoch at which the braneworld universe (B2) begins to accelerate. 
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FIG. 8. Confidence levels at 68.3% (light grey inner contour) 95.4% (medium grey contour) and 
99.73% (dark grey outer contour) in the n^^-^jn plane for the 'Disappearing Dark Energy' (DDE) 
model . The dashed lines show the current effective equation of state: wq = —0.48, —0.44, —0.4 
(left to right). 

in the BRANE2 model). 

Next we come to the third braneworld model, namely 'Disappearing Dark Energy' (DDE) 
for which the confidence levels in the Qm — plane are shown in figure 8. This case is 
different from the others in that the best fit value of flm tends to zero for this model. This 
behaviour is easily understood in view of the fact that the effective equation of state in this 
model (21) is always constrained to be >, — 0.5 and therefore cannot be very negative. 
Since the luminosity distance increases as flm decreases, a sufficiently large value of the 
luminosity distance (which agrees well with SNe observations) can only be achieved at the 
expense of having a small value of Qm in our braneworld model. Our results, shown in 
figure 8, are for the DDE model constrained by the prior Qm > 0.2. The best fit DDE model 
{Qjn — 0.2, — 0.8) has a chi-squared per degree of freedom of Xdoi — 1-20 which, although 
larger than the value for LCDM (xiof = 1-06 for = 0.3, x^^f = 1-08 for = 0.2), is 
much smaller than that for SCDM (Xdof — 1-75). 

18 




in 




0.5 1 1.5 

redshift z 



FIG. 9. The angular-size distance (in units of Hq ) in the braneworld models and in 
LCDM, SCDM is shown as a function of the redshift. The model parameters are: BRANEl 
{^m = 0.3, ni = 0.3, = 0), BRANE2 {Qru = 0.25, = 0.1, = 0), LCDM 
{rirn = 0.3, f^A = 0.7). All three models satisfy existing supernovae bounds. Also shown is the 
standard cold dark matter model (SCDM) with Qm = 1- 
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FIG. 10. The dimensionless quantity H{z)dA{z) which plays an important role in the Al- 
cock-Paczynski test is shown for braneworld models and for LCDM, SCDM as a function of the 
redshift. The model parameters are the same as in the previous figure. 
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FIG. 11. The quantity d\H~^{z) which plays a key role in the volume-redshift test is shown 
for braneworld models and for LCDM, SCDM as a function of the redshift. The model parameters 
are the same as in the previous figures. 
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FIG. 12. The dimensionless statefinder r = d /aH^ is shown for braneworld models and for 
LCDM. Also shown is the quintessence model with equation of state w = —2/3. The model 
parameters are the same as in the previous figure. 




redshift z 

FIG. 13. The dimensionless statefinder s = (r — l)/3(g — 1/2) is shown for braneworld models 
and for LCDM. Also shown is the quintessence model with equation of state w = —2/3. The model 
parameters are the same as in the previous figure. 
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One may try to constrain braneworld models further using some of the cosmological 
tests outhned earher in this paper. Through these tests it may be possible to distinguish 
between different braneworld models which agree well with current SNe data. To illustrate 
this we choose two complementary braneworld models: (a) BRANEl with = 0.3, Qi — 
0.3, Qa, = 0, and (b) BRANE2 with = 0.25, Q/ = 0.1, Qa, = 0, which satisfy both the 
high-^ SNe observations (at the la level) and the observational constraints on O^. In figure 
(9) we show theoretical results for the angular-size distance (Ia, which is determined from 
observations of the angle subtended by a 'standard ruler' of length £ at a redshift z, 
A9 = i/dA{z). Another useful quantity is H{z)dA{z), shown in figure (10), which plays a 
crucial role in the Alcock-Paczynski anisotropy test. We clearly see that the angular-size 
distance provides a better test with which braneworld models may be distinguished from 
each other and from rival models such as LCDM and SCDM. The Alcock-Paczynski test 
may not serve this purpose, since both Bl & B2 give results which are close to each other 
and to those for LCDM. The reason for this seeming degeneracy hes in the fact that, for Bl 
(B2), the Hubble parameter H{z) decreases (increases) relative to LCDM while the distance 
dA increases (decreases). Thus H{z) and show opposing tendencies for the braneworld 
models, which virtually neutralise each other in the combination H{z)dAiz). Exactly the 
opposite effect is achieved by the combination d\{z)H~^{z) which plays a key role in the 
volume-redshift test [19]. In this case the difference between Bl and B2 becomes more 
pronounced as demonstrated in figure 11. We therefore conclude that the volume-redshift 
test is probably a better means for differentiating between the two braneworld models Bl 
& B2 than the Alcock-Paczynski test. 

Interestingly the statefinder statistic (18) may also provide us with a useful means by 
which to discriminate between rival braneworld models and LCDM/quintessence, as demon- 
strated in figures (12) & (13). Other cosmological tests which are likely to be useful in testing 
braneworld models include gravitational lensing [24] and the cosmic microwave background 
[25]. We shall return to these important issues in a future work. 
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IV. DISCUSSION AND CONCLUSIONS 



This paper examines braneworld models of dark energy in the hght of recent supernova 
observations which indicate that the universe is accelerating. The braneworld models which 
we examine in this paper have several interesting properties which distinguish them both 
from the cosmological constant as well as from scalar field based 'tracker' models of dark 
energy. Like the latter, braneworld models presently accelerate, and possess a longer age 
than the standard cold dark matter model (SCDM). However in marked contrast to both 
LCDM and tracker models, the luminosity distance in one class of braneworld models, 
Bl, can be greater than the luminosity distance in LCDM (for identical values of VLm)'- 
di^iz) > d^^^^^^{z) > d^'^^^{z), where d^^{z) refers to the luminosity distance in 
spatially flat de Sitter space. In terms of the effective equation of state w, this is equivalent 
to the assertion that w < —1. This result is particularly surprising since matter in the 
braneworld model never violates the weak energy condition p+p > 0. A maximum likelihood 
analysis which compares braneworld model predictions with high redshift type la supernovae 
data, shows that Bl models provide good agreement with observations if Qrn ^ 0.3. These 
results broadly support the analysis of [23] in which 'phantom' dark energy models, having 
the property w — P/ p < —1, were compared against supernova observations. (It should be 
pointed out, however, that 'phantom models' invariably run into a physical singularity in 
the future when Pphantom — oo, such singularities are absent in the Bl model which remains 
well behaved at all future times.) 

The second braneworld model we consider (B2) has properties which complement those 
of Bl, since d^^^^^z) > d^^^^^'^{z) > d^^^^{z). This is equivalent to the assertion that 
—l<w<0. Results of a maximum hkelihood analysis show that B2 models are in excellent 
agreement with SNe data for smaller values of the density parameter fi^ <^ 0.25. 

Finally braneworld models also permit the dark energy to be a transient phenomenon. In 
models of this kind (called disappearing dark energy: DDE) the acceleration of the universe 
takes place during a transient regime separating past and future matter dominated epochs. 
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In these braneworld models, the universe does not possess an event horizon and so it may be 
possible to reconcile a universe which currently accelerates with the demands of string/M- 
theory. Comparison with Sne bounds shows that the Disappearing Dark Energy models 
marginally satisfy existing supernova data provided is sufficiently small: <, 0.23. 
For larger values of fi^, this class of models may be on the verge of being ruled out. 
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